Probing interactions in mesoscopic gold wires 
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We have measured in gold wires the energy exchange rate between quasiparticles, the phase 
coherence time of quasiparticles and the resistance vs. temperature, in order to probe the interac- 
tion processes which are relevant at low temperatures. We find that the energy exchange rate is 
higher than expected from the theory of electron-electron interactions, and that it has a different 
energy dependence. The dephasing time is constant at temperatures between 8 K and 0.5 K, and 
it increases below 0.5 K. The magnetoresistance is negative at large field scales, and the resistance 
decreases logarithmically with increasing temperatures, indicating the presence of magnetic impu- 
rities, probably Fe. Whereas resistivity and phase coherence measurements can be attributed to 
magnetic impurities, the question is raised whether these magnetic impurities could also mediate 
energy exchanges between quasiparticles. 
Several recent experiments have demonstrated that the low-energy properties of quasiparticles in metallic thin 
films are sample-dependent. On the one hand, the power-law increase of the phase-coherence time with decreasing 
temperature, predicted by the theory of electron-electron interactions in diffusive wires [jp, has been observed in 
several experiments J^|| . The energy exchange rates between electrons was found to be in agreement with this theory 
in experiments on silver wires [Q. On the other hand, the dephasing rate of quasiparticles was found to saturate at 
low temperature in a series of gold wires [^j , and the energy exchange rate between quasiparticles in copper wires to 
display an energy dependence different from the predicted one We present here measurements on gold samples, in 
which the energy exchange rates have the same energy dependence as was observed in copper and have a magnitude 
even higher. In order to investigate the origin of this effect, we have performed resistance measurements on samples 
fabricated similarly. The logarithmic dependence of the resistance, the negative magnetoresistance at large field and 
the temperature dependence of the phase coherence time, which is constant between 8 K and 0.5 K and increases at 
lower temperature, suggest the presence of magnetic impurities, which might mediate electron-electron interactions. 



I. ENERGY EXCHANGE RATES 
A. Measurement set-up 

In order to access the energy exchange rates among quasiparticles, we have measured the distribution function 
f(E) in a stationary out-of-equilibrium set-up described in Fig. |l|. We consider a mesoscopic metallic wire 
placed between two ideal reservoirs. In the absence of interactions, the population of quasiparticles at a given energy 
interpolates linearly between the distribution functions in the contacts, leading, if ksT <C eU, to a double-step- 
shaped distribution function. In the opposite "hot electron" regime, in which the typical interaction time is much 
shorter than the diffusion time Tr> — L 2 /D, equilibrium is reached locally at each position along the wire: the 
energy distribution function f(x,E) is a Fermi function (see dotted curves in Fig. Q). If heat is only carried out 
by electrons, the local temperature is T e g(x) = yjT 2 + x (1 — x) U 2 /£, where £ = (^f-) is the Lorenz number 
10,^). Our experiments focus on the intermediate regime, in which interactions lead to a significant redistribution 
of the energy among quasiparticles, but not to a complete thermalization. The distribution function is obtained 
from the differential conductance dl/dV (V) of a tunnel junction between the wire and superconducting electrodes 
HI). When the temperature of the superconductor lies well below the superconducting transition temperature, and 
if the density of states in the normal electrode wire is taken as energy independent on the probed energy range, the 
differential conductance of the junction is simply proportional to the convolution product of the density of states in 
the superconductor ns(E) — \E\ /V E 2 — A 2 , with A the gap of the superconductor, and of the energy derivative of 
the distribution function in the wire, 3 ^g^" E ' ) : 
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with R t the tunnel resistance of the junction. 
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FIG. 1. Experimental layout: a metallic wire of length L is connected to large reservoir electrodes, biased at potentials and 
U. In absence of interaction, the distribution function at a distance X = xL from the grounded electrode has an intermediate 
step / (E) = 1 — x for energies between —eU and (solid curves) (we assume U > 0). When interactions are strong enough to 
thermalize electrons, the distribution function is a Fermi function (dotted curves). In the experiment, the distribution function 
is obtained from the differential conductance dI/dV(V) of the tunnel junction formed by the wire and a superconducting 
electrode placed underneath. 



B. Samples 

The samples were fabricated by deposition with an electron-gun, at several angles, through a PMMA suspended 
mask patterned using e-beam lithography. The substrate was thermally oxidized silicon, as in our experiments on 
Cu and Ag. We fabricated on the same chip two Au wires: wire #1, with length L — 1.55 /im, and a single probe 
junction placed at x = 0.7; and wire #2, with length L = 5 /im, and two probe junctions, at x = 0.25 and x = 0.5. 

We first deposited a 25 nm-thick aluminum film, which was then oxidized. This layer defines the superconducting 
probe electrodes. The wires and the pads were obtained by the subsequent evaporation from a 99.99% purity gold 
target, at a pressure of 10 -6 mb, at 1 nm/s. The thickness and width of the wires are 45 nm and 165 nm. The 
electrodes at the ends of the wires are 500 nm-thick pads, with an area of about 1 mm 2 . From the low temperature 
wire resistances R\ = 5.39 VL and R2 = 16.9 f2, we deduce from Einstein's relation, assuming rectangular cross-sections, 
the diffusion constant D — 0.013 m 2 /s and the diffusion times tdi = 0.18 ns and td2 — 1.8 ns. The samples were 
mounted in a copper box thermally anchored to the mixing chamber of a dilution refrigerator. Electrical connections 
were made through filtered coaxial lines |10|. 
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FIG. 2. Micrometric scale (left) and large scale (right) shape of the sample. The wires are placed between pairs of large 
contact pads (wire #1 between A and B, wire #2 between A and C). Superconducting probe junctions are placed on wire #1 
at x — 0.7, on wire #2 at x — 0.25 and x — 0.5. The shape of the pads A, B, C has been designed to optimize the cooling of 
the reservoirs. 

Measurements proceed as follows. In a first step, the voltage U is set to zero. From the comparison of the measured 
differential conductance of the tunnel junctions with the calculated convolution product of the BCS density of states 
and of the derivative of a Fermi function at the temperature of the mixing chamber (see Eq. (1)), we deduce the gap 
of the superconductor (A w 0.2 mV) and the tunnel resistances |flj (Rt = 57 kfl (wire #1), 28 kf2 and 30 kf2 (wire 
#2)). In a second step, the voltage U is set to a dc value. From the differential conductance of the tunnel junctions, 
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we deduce by numerical deconvolution the distribution functions in the wire at the position of the junctions, using 
Eq. (1) and the values of the gap and tunnel resistance determined in the first step. 

C. Distribution functions 

We show in Fig. || the distribution functions measured with the three junctions at U = 0.1, 0.2, 0.3 and 0.4 mV 
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FIG. 3. Measured distribution functions f(E) on wire #1 (left panel) and wire #2 (central and right panel) for 
U — 0.1, 0.2, 0.3 and 0.4 mV, plotted as a function of E (top) and E/eU (bottom). The dotted line in the central panel 
is the prediction in the hot electron regime for U = 0.4 mV. 

As expected from the difference in the diffusion times, the distribution functions are much more rounded in wire 
#2 than in wire #1. The distribution function at the middle of wire #2 is close to the hot electron prediction, shown 
with a dotted line in the central panel for U = 0.4 mV. However, the distribution functions measured on the same 
wire close to the left electrode (right panel) are clearly different from Fermi functions: they display a strong slope 
near E = 0, reminiscent of the Fermi step at this energy in the closest pad (A) which was grounded. 

The larger the voltage U, the wider the interval of energy in which / varies from 1 to 0, as expected in both limiting 
regimes (see Fig. |l|). In order to remove this dependence, we have replotted the same data in the bottom of Fig. ||, 
but with the reduced energy E/eU on the horizontal axis. Remarkably, all scaled curves superimpose. This scaling 
property, which was also found in all the experiments on copper wires, is not generic: in experiments on silver wires, 
the slope of f(E) at E/eU — —0.5 was found to increase with U Q. 
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D. Energy exchange kernel 

The distribution function can be calculated by solving the stationary Boltzmann equation in the diffusive regime 
©|: 

^^qI^ + Z con E, {/}) - 2££ (x, E, {/}) = (2) 

where ZJ." U (x,E, {/}) and 2"°")* (x, E,{f}) are the rates at which quasiparticles are scattered in and out of a state 
at energy E by inelastic processes. The observation of the hot-electron regime in the middle of wire #2, with a 
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temperature close to the expected one, indicates that the energy is mainly redistributed among the quasiparticles. In 
particular, phonon emission can be neglected. Assuming that the dominant inelastic process is a two-quasiparticle 
interaction which is local on the scale of variations of the distribution function, 

T» n (x, E,{f}) = J dsdE'K(s)r E+s (l -f%)f%, (l-fp_ e ) (3) 

where the shorthand /Ji stands for / (x, E) . The out-collision term 1°^' has a similar form. The kernel function 
K (e) is proportional to the averaged squared interaction between two quasiparticles exchanging an energy e. We have 
neglected the possible dependence of K(e) on the energies of the initial and final states and on the position along the 
wire. The theory of electron-electron interactions in diffusive conductors in the ID regime Q predicts K (e) cx e -3 / 2 , 
a regime which was observed in silver wires [Q. In gold and copper wires, the scaling property implies, by a simple 
change of variables in Eq. (3), that U 2 K (e) is a function of e/eU only |Q. If futhermore K (e) does not depend on U, 
one obtains K (e) = j/e 2 , with 7 a typical interaction rate. We nevertheless tried to fit our data at U = 0.1 mV with 
K (e) = k/e 3 / 2 , and obtained k = 50 ns _1 meV -1 ^ 2 , which is three orders of magnitude larger than the theoretical 
prediction K thy = {fc'^vSit y' D /2)~ 1 sa 0.06 ns _1 meV -1 ^ 2 , with v the density of states in gold and S the cross- 
section of the wire Jl5|,p[. Moreover, the shape of the distribution functions on wire #2 is not properly reproduced 
and the calculated curve for wire #1 at U = 0.4 mV with the same parameter is significantly more rounded than the 
experimental data. 



E. Fits 



Figure |J shows the best fit of the data with K (e) = 7/e 2 , obtained with 7 = 8.9 ns 







1 1 1 ! 1 


1 1 1 1 1 


1 1 1 1 1 


0.4 mV 


"~~^X\ 0.4 mV 


0.4 mV 


0.1 mVV^^ 


0.1 mv\\ 


0.1 mv\\ 




wire #1 ^\ V___ 


wire #2 \ 


wire #2 \ x,^ 


x-0.7 \ 


x-0.5 \. 


x=0.25 


i.i. 


i.i. 


i.i. 


-1 


-1 


-1 


E/eU 


E/eU 


E/eU 



FIG. 4. Symbols: measured distribution functions with U = 0.1 mV and (shifted vertically) U = 0.4 mV. Solid lines: calcu- 
lated distribution functions with the phenomenological kernel K (e) = 7/e 2 . Reservoirs are here assumed to be at T — 33 mK. 



By construction, the scaling property of the data is exactly reproduced. Whereas the shape of the distribution 
functions on wire #2 is well accounted for, the experimental data on wire #1 are more rounded than we calculate. 
The opposite discrepancy was observed in experiments on 5 /im— long copper wires fl6|| . We attribute the strong 
rounding in the distribution functions on wire $d to the large heating of the reservoirs associated with the low 
resistance of wire #1: at a given value of voltage U, the heat power P = U 2 / R is the highest in the less-resistive 
wires. The solution of the heat equation in the contacts, assumed to have the shape of an angular sector, with angle 8, 
and neglecting phonon emission, gives the reservoir temperature at the end of the wire ]17j : T = \/Tq + b 2 U 2 . Here, T 

is the temperature of the quasiparticles at a large distance r max from the contact to the wire, and b = J — ^- In , 
with Rq 8 pa 0.05 fl the sheet resistance of the reservoir, estimated from the resistivity of the wires and the thickness 
ratio of the wires and of the reservoirs, r max « 1 mm a typical equilibration length between electrons and phonons 
p8| , and r m i n the smallest radius for the radial approximation to be valid: r m i n « w. In the experiment, the left 
contact (labelled A in Fig. ||), opens with an angle 6l ~ 2.44 rad, whereas the two other contacts (B and C in Fig. ||) 
have a smaller angle: Or ~ 0.96 rad. We have fitted the distribution functions separately from U = 0.005 mV to 
U = 0.5 mV, with the temperatures of the reservoirs as fit parameters. The results are shown in Fig. H. 
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FIG. 5. Left panel : distribution functions on wire #1 calculated with K (e) = 7/e 2 at U = 0.1 mV and 0.4 mV, the reservoirs 
temperatures being taken as fit parameters (solid lines). Symbols are data points. Right panel: symbols: temperatures of the 
right (Tr) and left (Tl) reservoirs obtained from the fits. Solid lines: reservoir heating model with b as a free parameter. 



At voltages below 20 /j,V, the fit temperature is nearly constant and identical in both reservoirs: T = 43 mK, in 
reasonable agreement with the temperature T = 33 mK indicated by the thermometer on the mixing chamber. At such 
low voltages, scaling is not obeyed since temperature produces a significant rounding of the steps and eU /ksT varies 
with U. From the fit of the data on Fig. ||, we deduce 6l = 0.50 mK//xV for the left contact, and bn = 0.80 mK//xV 
for the right contact. Theory gives b^jb^ — \/&l/9r, = 1.60, in excellent agreement with the ratio found from the 
fits. Our estimated values of 6l and bn are about a factor of 2 larger than those obtained from Fig. [|, possibly 
because the sheet resistance of the reservoirs is less than our estimate based on the electronic mean free path in the 
wire. At voltages larger than 50 /iV, the temperature of the reservoirs is proportional to U, leading to a scaling in the 
rounding of the steps. In contrast, we want to stress that the slope of the distribution function near E/eU = —0.5 
cannot be explained by heating alone. The fits of the distribution functions in this regime, taking into account both 
K (e) = 7/e 2 (with 7 = 8.9 ns" 1 ) and heating is shown on the left panel of Fig. ||. A consequence of reservoir heating 
is that the distribution functions are rounded at the scale of ^f-bU ss U/15, and the sharp features that were observed 
on copper samples cannot not be resolved [lq| . 



F. Two- level systems 



The distribution functions can also be well accounted for by a model in which quasiparticles are in local equilibrium 
with two-level systems distributed uniformly along the wire. The relevance of TLS on phase relaxation has recently 
been suggested by several authors ]l9|-|2l[|. Agreement is found if one assumes, in a simple model (described in 
that the quasiparticles are weakly coupled to the TLS with a density inversely proportional to the spacing e between 
the two levels. Such a density is obtained if the two-level systems are the two lowest energy levels in symmetric 
double- well potentials, and if the distribution of the barrier heights is uniform. In glasses the distribution of potential 
well asymmetries is usually taken as white p^ ], but one might argue that in metals, symmetric double- wells could 
result from crystalline symmetries p3[ . We note that a similar assumption is made in calculations based on a strong 
coupling to TLS 0. 



II. RESISTANCE MEASUREMENTS 



Complementary information on interactions was obtained from resistance measurements. We have fabricated long 
gold wires in the same deposition machine as for the energy relaxation experiment. To enhance adhesion of the Au film 
on the substrate, we used two different methods. For a first sample (Aul) we evaporated first 1 nm of aluminum, 
and oxidized it. On a second sample (Au2), the surface of the sample was ion-milled just before gold deposition. A 
more complete set of data was taken on sample Au2, and we report only here the results on this sample. Its length, 
width and thickness are L = 271 /im, w = 115 nm, and t = 45 nm, respectively. The low temperature resistance was 
1125 ft. Assuming Einstein's relation and a rectangular cross-section, we deduce the diffusion constant D = 0.016 m 2 /s 
. In another fabrication run conducted at Michigan State University, we have Joule evaporated another sample, 
called AuMSU, with very pure gold (99.9999%, i.e. 1 ppm of impurities), with L — 176 /xm, w = 80 nm, t = 45 nm 
and D = 0.016 m 2 /s. All the samples were measured in the same top-loading dilution refrigerator. 
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A. Resistance vs. temperature 



We show in Fig. || the temperature dependence of the resistance for samples Au2 and AuMSU. The contribution 
of weak localization for Au2, smaller than 10 -4 (see below), can be neglected. For AuMSU, this contribution, up to 
3x 10~ 3 in SR/R, has a well understood temperature dependence (see below), and was substracted. The contribution of 

electron-electron interactions to the resistivity, ~ 3.13-^^ = -j= p5[ ], with Rk = an d Lj- = yQ^", plotted 

as dotted lines in Fig. ^, accounts well for the data of AuMSU pf|, where the fit parameter a = 2.7 x 10~ 3 K^ 1 / 2 is 
very close to the calculated value a thy = 2.5 x 1CP 3 K^ 1 / 2 . In contrast, the variations observed in Au2 are stronger 
and have a different temperature dependence. 
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FIG. 6. Relative variations of the resistance of sample Au2 and AuMSU as a function of temperature. The dotted line is the 
prediction of the theory of electron-electron interactions; the solid line on the left panel is a fit with an additional logarithmic 
term. 



A good fit of the data can be obtained with: ^ = — (3\ogT, with (3 = 1.23 x 10~ 3 . The logarithmic term 
is characteristic of the Kondo effect 1 27 , and it was in particular studied in great detail in dilute alloys of Fe in Au 
p8| pl| , where (3 was found to be proportional to the Fe concentration and dependent on the cross section of the wire 
as well as on disorder. Assuming that the logarithmic dependence of the resistance of Au2 is due to Fe impurities, we 
get, by comparison with a sample with similar width, thickness and diffusion constant (AuFe2 in ]3l[|), an impurity 
concentration c w 55 ppm, which is compatible with the purity of the gold source used for evaporation (99.99%, i.e. 
100 ppm of impurities) |B3]. 



B. Resistance vs. magnetic field - Phase coherence time 

We present in the left panel of Fig. |7] the magnetoresistance of Au2 and AuMSU. Around zero magnetic field, 
the resistance presents a dip, as expected from weak antilocalisation [jlj. A large-scale negative magnetoresistance is 
found on Au2, which is similar to measurements performed on dilute alloys of Fe in Au pS|J3l]]. It can be attributed 



to the freezing of the magnetic moment of the impurities with the magnetic field |28|. In order to extract the phase 
coherence time from the magnetoresistance of Au2, we subtracted out the large scale magnetoresistance, which was, 
at each temperature, fitted with AR(B)/R = -r(T)(yJ B 2 + B 2 - B ). The remaining magnetoresistance was then 
fitted with the predictions of the weak localization theory [|||- In the right panel of Fig. 0, we show the phase 
coherence time r v of samples Aul, Au2 and AuMSU as a function of temperature. Upon cooling, the phase coherence 
time of Aul and Au2 remains unchanged from 8 K to 0.5 K: t v ps 10 ps, then it increases roughly as 1/T. A simple 
interpretation of the very low value of t v and of the desaturation below 0.5 K is found by comparing with similar 
measurements performed on gold wires in which a small amount of iron was purposely introduced ]2^,^9|,[5],[3l] . The 
same behavior of t v (T) was found in these experiments, and it was interpreted as an effect of spin-flip scattering by 
the impurities, the rate of which presents a maximum at the Kondo temperature of Fe in Au (Tr- 0.3 — 1 K). The 
value r° of t v near the Kondo temperature was found to be roughly given by r° « 0.25 ns/c, with the concentration of 
impurities c expressed in ppm (parts per million) [5|j3l[|. In this interpretation, « 10 ps is obtained for c ~ 25 ppm. 
This provides us with another estimation of c which is of the same order of magnitude as what we deduced from 
the resistance vs. temperature measurements. Sample AuMSU does not show any magnetoresistance at large scale, 
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and the phase coherence time presents no saturation down to 44 mK (see Fig. [?]). A good fit of the data could be 
obtained with the theoretical dependence t^ 1 = AT 2 / 3 + BT 3 , with A = 0.9 ns^K^ 2 / 3 and B = 68 /xs _1 K- 3 . 
The theoretical value from the theory of electron-electron interactions is A = 0.5 ns _1 K~ 2 / 3 [f25f . 
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FIG. 7. Left panel: Symbols: magnetoresistance of sample Au2 (vertical shifts between the curves taken at different tem- 
peratures are arbitrary) and of sample AuMSU (inset; note the difference in the scales). Solid curves: fit with the sum of 
the weak localisation expressions and (for Au2 only) of a phenomenological large scale contribution —r(T)^/Bg + B 2 . Right 
panel: phase coherence time as a function of temperature, for samples Auf , Au2 and AuMSU. Solid line: fit of AuMSU with 
r" 1 = AT 2/3 + BT 3 . For reference, we have also plotted the data of two samples of Ref. [[| similar to AuMSU (MJW Au-3&6) 

To allow for a comparison with the results of Mohanty et al, we have also plotted on the same hgure the phase 
coherence time of samples Au-3 and Au-6 in 0, which have similar size and diffusion constants as AuMSU (Au-3: 
W = fOO nm, t = 35 nm and D = 0.011 m 2 /s; Au-6: w = 180 nm, t = 40 nm and D = 0.016 m 2 /s). It is clearly seen 
that our data do not display the saturation found in [^| . We had already found no saturation in measurements on silver 
wires ||, which were made from very pure (99.9999%) Ag. We suspect that energy exchanges in samples fabricated 
with very pure Au would be consistent with the theory of Ref. Q , as in very pure Ag , but this experiment remains 
to be performed. It is noteworthy that the low-temperature saturation of reported in Cu wires |3j was identically 
found in other samples obtained from 99.9999% Cu, supporting the hypothesis that surface oxide could play a major 
role in Cu @. 

III. CAN MAGNETIC IMPURITIES MEDIATE ELECTRON-ELECTRON INTERACTIONS? 

The data presented in section 2 provide strong evidence that some of the gold wires we have studied contain several 
tens of ppm of iron. The question arises as to whether the energy exchange measurements can also be explained 
by the presence of these magnetic impurities. Kaminski and Glazman |3^| recently pointed out that K (e) oc e~ 2 is 
obtained in second order perturbation theory from the interaction of two quasiparticles with a magnetic impurity. The 
effective matrix element M e ff of the second order process is proportional to J 2 /(E V — Ei), where J is the coupling 
parameter between the quasiparticles and the magnetic impurity and E v — E% is the energy difference between the 
intermediate (virtual) state and the initial state. In the intermediate state, one quasiparticle is promoted to an energy 
state higher by e than the initial state, whereas the magnetic impurity has been reversed. One thus obtains K(e) oc 
c x |Af e yy| 2 = cJ 4 /e 2 , which has the energy dependence observed in the experiment, but the bare coupling parameter 
J w l/v\og(Ep/kBTx) is too small. If one assumes that we have Fe impurities in Au, kgTx ~ 50 /LteV is of the order 
of the energies probed in the energy exchange experiment, and the coupling constant is expected to be renormalised by 
the Kondo effect. This effect has been treated in [Q, and the proposed expression K(e) = ^j^S(S + 1) log _4 (|^-)p-, 
has now the right order of magnitude but too strong a U -dependence. Further work is clearly needed to clarify this 
issue. 
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